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Abstract 

We investigated numerically dyon-like solutions of the SU(2) Einstein- Yang-Mills system on a cylin- 
drically symmetric space time with a cosmological constant. We find a new kind of behaviour not found 
in the spherically symmetric models. For positive values of A we have an oscillatory behaviour of the 
magnetic component of the YM field around the r-axis, so there is an arbitrary number of nodes. 
For increasing positive A, the frequency increases also and the solution breaks down at finite r = r^, 
indicating a singularity. The electric component, however, approaches a constant value. After further 
increasing A, this global behaviour repeats itself at a larger r while the former singular behaviour 
disappears. This behaviour depends critically on the position of the core of the string. For increasing 
negative A, the oscillatory behaviour disappears and the magnetic and electric components behave like 
the scalar and gauge field in the Abelian cosmic string model. 

1 Introduction 

It is known for decades that the static Einstein- Yang-MiUs theory for the gauge group SU(2) admits 
regular solutions. The first example was that of Bartnik and McKinnon (BK) The basic idee behind 
this phenomenon is the balance between the attractive gravity and the repulsive Yang-Mills field. When 
a cosmological constant is incorporated in the model, the behaviour of the particle-like solution can be 
significant different compared to the original BK solution [21 El E El- For a nice overview, see 

When the cosmological constant is negative, there are two types of generic solutions, the singular ones 
and the asymptotical anti-de Sitter (AdS) solutions. There is also a change of the topology of the moduli 
space of globally regular ('soliton') and black hole solutions when going from negative to positive values 
of A. This is called a 'fractal' structure 3^. Moreover, it is conjectured that negative A solutions are the 
only stable EYM solutions and the magnetic component of the YM field need not be nodeless in order 
to be stable. In these models, one usually takes the pure magnetic ansatz. Further, one can investigate 
the dyon- type solutions, where a non-zero electric component of the YM field is incorporated I2j. It turns 
out that dyon solution can exist in these models with negative A, however under quite special asymptotic 
assumptions. In the present paper we investigate the EYM system with positive and negative cosmological 
constant on a cylindrical space time, without any assumption on the asymptotic behaviour and try to find 
a connection between the Abelian sting-like solution found on the same space time. 
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2 The field equations of the Einstein Yang-Mills system 

Consider the Lagrangian of the SU(2) EYM system 



S 



d X\J—g 



with the YM field strength 



(1) 



(2) 



g the gauge couphng constant, G Newton's constant, the gauge potential, and TZ the curvature scalar. 
The field equations then become 

G^i/ = SttGTi^u - Ag^i., (3) 

(4) 

(5) 



with T the energy-momentum tensor 



'77a ^af3a 



and V the gauge-covariant derivative. 



Let us consider the stationairy exterior line element 

ds^ = C{-dt^ + dz^) + dr^ + B^r^difi^, 
where B and C arc functions of r. For the YM-field we will take [3IH| 

A^dxf" = ($if^ + <p2Tz)dt + AiT^dr + A2t^dz + {Wifr + {W2 - l)f^)dip 
with $i, Ai and Wi functions of r, and ti the orthonormal axial generators of the SU(2) 



(6) 
(7) 
(8) 



r 



{cos{mp), sin{n(p), 0), 



= (0,0,1), 



(sin{n(l)), — cos(ri(/)), 0), 



(9) 



with n the magnetic charge. For the special case, Ai — A2 — 0, $2 = Wi — and n=l, we obtain a 
consistent set of equations: 



C" 



1 



l\2 



AC ^ 



AC, 



q,^W^ + B^r\^'Y - 2C{W') 



CW 
Br 



Br ' 
Br ' 



where k = The conservation law V ^T^'^ is also fulfilled. Further, there is a scale invariance 



C ^ 7C, $ ^ ^7*' B^ PB, W ^ f3W, 



(10) 



(11) 



with P and 7 some constant. 
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3 Numerical solution 



The set of equations are easily solved numerically. We used two different Fortran ODE-solvers to check 
the reliability of the numerical solution and the MAPLE-program as final check. For all the routine-codes 
we find the same solution. The absolute and relative tolerance was I.IO"'^^ For some specific values of 
"shooting" parameters , we find a remarkable solution critically dependent on the cosmological constant. 
Let us consider near the core of the string ro : 



W^(ro) = l-K + ..., $(ro)=aro + 
Then we obtain for the initial values of the metric components: 



0{rl), C{ro) = 1 - x/31n(3 - f .5^3) + 6)+0(r^) 



(12) 



(13) 



In the spherical symmetric situation ( |21) the different solutions are classified by some parameters de- 
termined by the condition of the field variables at r — > oo. Here we are mainly interested in the critical 
situation for some initial values at rg The energy density becomes 



l\2 



(14) 



We observe that for the chosen initial values, Ttt is finite for r tq. 

In figure 1 and 2 we plotted a typical solution for positive and negative values of A, and for a — 
1,6 = 0.35 and for tq = 0.01. For positive values of A we observe an oscillatory behaviour of the magnetic 
component around zero, so there are an arbitrary number of zero's. For increasing positive A, the frequency 
increases also and the solution breaks down at finite r = Vh where C approaches zero. This means a 
cosmological horizon. However, for some value of A ( for example A around 0.004) the horizon is pushed 
to larger values of r, while the former horizon is resolved. This happens again for larger positive values of 
A. In figure 3 we plotted the horizon r/j as function of A. For increasing negative A we observe that the 
oscillation of W slows down and disappears at a critical A. The solution behaves like a string solution. 
This fact was also observed in the radiative counterpart model [S] . For smaller values of ro this pushing 
foreward of r/j dissapears. In figure 4 we plotted Vh versus A for ro = 10^^. 
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Figure 1 
respectively. 



Plot of B{t),C(t),W{t), and $(r) for negative A. We took for A -0.25, -0.15, and -0.01 
= 1 and 6 = 1, c = 0.35.ro = 0.01. 
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Figure 2. Plot of B{r),C{r), W{r), $(r) for positive A and initial values a = l,b = 0.35. We took A 0, 
0.0005, 0.001, 0.0025, 0.0034, 0.00425, 0.0075, 0.0125, 0.014, and 0.015 respectively, k = l.ro = 0.01. 
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Figure 3. Plot of the horizon versus A, for ro = 0.01. 
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Figure 4. Plot of the horizon r/j versus A, for ro = 10 



4 Conclusions 

This paper is a first attempt to investigate numerically the dyon solutions of the EYM model with a 
cosmological constant on a cylindrically symmetric space time. In the numerical code, we start with 
regular initial values of the field variables at a small value of the core radius of the string and varied A 
continuously positive as well as negative. 

In contrast with the spherical symmetric EYM solutions with negative cosmological constant, our 
numerical solutions seems to be regular every where without any a priori assumption on the asymptotic 
behaviour. The solution behaves like a cosmic string. However further investigation is needed. For 
positive cosmological constant, the numerical solutions are quite different from the spherically symmetric 
counterpart models. First of all, the rapid oscillatory behaviour of W is not found in earlier studies. The 
classification on the number of nodes of W makes no sense in our solution. Secondly, The position of Vh 
moves towards the origin for increasing A, but shows a periodic behaviour when A is further increased. 
When the core radius is decreased, this periodic behaviour dissapears. The solution found here needs 
further investigation, specially for different initial value. Further, the charges should be incorporated in 
the code, for example, the electric 

QE = ^{j r^Wr^dr -r j $'r^d^) (15) 
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